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' A model of directed percolation processes with colors and flavors that is equivalent to a population 

model with many species near their extinction thresholds is presented. We use renormalized field 
theory and demonstrate that all renormalizations needed for the calculation of the universal scaling 
behavior near the multicritical point can be gained from the one-species Gribov process (Reggeon 
field theory). In addition this universal model shows an instability that generically leads to a total 
asymmetry between each pair of species of a cooperative society, and finally to unidirectionality 
QO ' °f the interspecies couplings. It is shown that in general the universal multicritical properties of 

unidirectionally coupled directed percolation processes with linear coupling can also be described 
by the model. Consequently the crossover exponent describing the scaling of the linear coupling 
parameters is given by $ = 1 to all orders of the perturbation expansion. As an example of 
^ ■ unidirectionally coupled directed percolation, we discuss the population dynamics of the tournaments 

of three colors. 
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I , I. INTRODUCTION 

Nonequilibrium processes, their stationary states and their phase transitions have been of considerable interest in 
natural science as well as in medicine and sociology for many years. Here we are interested in processes that can be 
modelled by growth and decay of populations with spatially local interaction rules. The transition between survival 
and extinction of a population is a nonequilibrium continuous phase transition phenomenon and is characterized 
by universal scaling laws. It is well known that also for systems far from equilibrium the concept of universality 
' classes with respect to their critical properties in the vicinity of a continuous phase transition is applicable. For 
the description of transitions in systems that show active and absorbing inactive states, percolation models play an 
outstanding role. Some years ago it was conjectured [QJ|] that Markovian growth models with one-component order 
parameters displaying a transition into an absorbing state in the absence of any special conservation law generically 
belong to the universality class of directed percolation (DP). Besides DP |^-|5| this universality class includes e.g. 
' Reggeon field theory (RFT) [|-||, the contact process P-[yJ , certain cellular automata |12j and some catalysis models 
-j-^ . [PfB (for a recent review of DP processes see |ll| ) . 

Despite the fact that a large variety of different models belong to the DP universality class, there is still no experiment 
where the critical behavior of DP was seen [ [l6| . In a recent paper |ljj , Hinrichsen compares suggested experiments 
and discusses possible reasons why the observation of DP critical exponents is obscured or even impossible. One of 
these reasons might be that the basic feature of the DP class, the existence of an absorbing state, is quite difficult to 
realize in nature. Small fluctuations will always affect this state and may be strong enough to soften the transition 
like a small particle source, which works as an external field Jl8| . Another reason might be the influence of spatial 
quenched disorder which is abundant in reality. We have shown ||l9|| that in contrast to equilibrium systems, the 
critical scaling properties of DP processes are not only altered by frozen randomness, but fully destroyed. 

For the analytic description of universal behavior near a critical nonequilibrium transition, it is often useful to model 
the universality class by mesoscopic stochastic processes involving the order parameter and other relevant fields. In 
03 1 case of the DP class a representation by the Langcvin equation for the time-development of the particle density, the 
Gribov process (the stochastic version of the so-called Schlogl model |2(|), is appropriate. The name Gribov process 
was coined by Grassberger who showed that RFT is a Markov process in disguise rather than a quantum theory Q . On 
the level of a formulation of stochastic processes by means of path integrals, there is superficially no difference between 
RFT and the Gribov process. However RFT uses creation and annihilation operators for particles as the principal 
fields in contrast to the particle density and its conjugate response field in the Gribov process. Microscopically the 
RFT-description of DP starts with special reactions between diffusing individuals on a given d-dimensional lattice 
such as birth: X — > X + X, competition: X + X — > X, and death: X — > 0. These reactions are represented by a 
master equation that is mapped onto a second-quantized bosonic operator representation, which is in turn mapped 
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onto a bosonic field theory using the continuum limit |2l[] ■ At the critical point the rates of birth and death have to 
balance to yield a vanishing overall production of individuals. But this condition leads to strong local correlations 
and the microstates of the system consist typically of clusters of individuals embedded in the vacuum 0). Thus, a 
fluctuating density description for the Gribov process is appropriate on a mesoscopic level. Note that the replacement 
of the spontaneous single particle death reaction by a two particle reaction X + X — > leads to strong anticorrelations 
because now the birth rate itself has to vanish at the critical point. In this case, the microstates consist of separated 
lonely wanderers that only sporadically interact and no clustering of individuals sets in. For such problems of branching 
and annihilating random walks one is forced to use the creation- annihilation operator formulation [p2[ and a naive 
stochastic density description for the fundamental field would lead to wrong results. As a general rule: one can show 
that branching processes lead to positive correlations and annihilation processes to anticorrelations with the exception 
of spontaneous decay, which does not generate any correlations. Thus, nonvanishing branching and spontaneous decay 
are needed to yield positive correlations together with the possibility of a vanishing overall production of particles at 
the critical point. In such cases a mesoscopic density description is correct. 

Instead of considering only one species of particles as is usually done for processes belonging to the DP class, it is of 
interest to introduce processes with several interacting species as for instance in mathematical biology [2^j2^| , which 
are also of relevance for a special model of surface growth [^5[ . It is the purpose of this paper to describe the detailed 
field-theoretic investigation based on renormalization group methods of such colored and flavored directed percolation 
processes (MDP) realized by the Gribov process for several species. We group different colored species in the same 
flavor class, if they have equal transport properties. 

In the next chapter we introduce the model and define its renormalization and the one-loop calculation in the 
third chapter. In chapter IV, we present the general renormalization group analysis and find the asymptotic scaling 
behavior to one-loop order. In chapter V, we include the two-loop results of the appendices and show the crossover 
to unidirectionality of the couplings between the different species. In chapter VI, we present our considerations on 
symmetries and the general fixed point properties of the model. We show that the permutation symmetry of a 
multicolored process is spontaneously broken. A brief account of this work has been presented in (26) . In addition, 
in chapter VII, we will show that the universal multicritical features of a recently introduced model of unidirectional 
coupled directed percolation processes (UCDP) by Tauber et al. ^7j, which contains an additional linear coupling 
between the species, is completely described by the MDP class. In the last chapter we discuss the results and give an 
application to biomathematics. Three appendices present technical details, e.g. the e-expansion of the DP-exponents, 
known for a long time, but as yet unpublished. 

II. THE MODEL 

The mesoscopic description of the dynamics of physical systems is based on a correct choice of the complete 
set of fundamental slowly-developing fields. In general these are the order parameter densities and the densities of 
conserved quantities. The multispecies processes under consideration are completely described by the particle densities 
n (x, t) = (n\ (x, t) , n.2 (x, t) , . . .) of the percolating colored and flavored individuals. We assume that there does not 
exist any conservation law. 

Next one has to find out the general form of the stochastic equations of motions of the fundamental fields as 
timelocal (Markovian) Langevin equations. These Langevin equations have to respect symmetries and general prin- 
ciples characterizing the universality class under consideration. The MDP-class is characterized by the following four 
principles: 

1. Errorfree self-reproduction ("birth") and spontaneous annihilation ("death") of individuals. The rates for birth 
and death may be different for each color. 

2. Interaction between the individuals ("competition", "saturation") with color-dependent couplings. 

3. Diffusion ("motion", "spreading") of the individuals in a <i-dimensional space with flavor-depending transport 
coefficients. 

4. The states with at least one extinct color are absorbing. 

In the language of chemistry, the MDP may be realized microscopically by an autocatalytic reaction scheme of the 
form X a «-> 2X a , X a — > 0, X a + Xp — » kX a + IX p, where the last reaction subsumes the interactions of individuals 
with colors a, f3, and k,l may be the integers 0, 1. A description in terms of particle densities typically arises from 
a coarse-graining procedure where a large number of microscopic degrees of freedom are averaged out. The influence 
of these is simply modelled by Gaussian noise-terms in the Langevin equation which however have to respect the 
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absorbing state condition. The stochastic reaction-diffusion equations for the particle densities in accordance with 
the four principles given above are of the form 



d t n a (x,t) = A Q V 2 n Q (x, t) + R a (n(x, t))n a (x,t) +Ca(x, t) , (1) 

where the first term on the right hand side models the (diffusive) motion, and the R a are the overall reproduction 
rates of the particles with color a. These deterministic terms are constructed proportional to n a in order to ensure 
the existence of an absorbing state for each species. Near the absorbing transition the particle densities n are small 
quantities. Expanding the rates R a in powers of n results in 

R a (n) = -A Q \r a + - H/H J. (2) 

The Gaussian noises £ Q (r, f) must also respect the absorbing state condition, whence 

(Ca(x,t)C /9 (x',0) = 2D aP (n(x,t))6{x - x.')5(t - t') + ■ ■ ■ 

= \ a g a 8 a ,pn a (x.,t)d(x.--x.')5(t-t') -\ . (3) 

Subleading terms in the expansions (^,[|) as well as additional terms with derivatives of the spatial 5-function in the 
first line of Eq. (^) are not displayed. It can be shown that they are irrelevant in the renormalization group sense 
as long as the stability condition „ \ a 9ai3na n i3 > f° r a U n Q > is fulfilled. The breakdown of this condition 
signals the occurrence of a discontinuous transition to compact growth and the appearance of tricritical phenomena 
at the border between first and second order transitions. Such a behavior is expected in microscopic models based on 
more complicated particle reactions p£| . The "temperature" variables r Q measure the difference of the rates of death 
and birth of the color a. Thus the temperatures may be positive ore negative. We are interested in the case where 
all T a « (up to fluctuation corrections) which defines the multicritical region. Under these conditions all the species 
live on the border of extinction. 

The next step is a mean field investigation of homogeneous steady state solutions of the equations of motion. 
Neglecting all fluctuations in Eq. ([!]) and using the expansion (g) up to first order in n, we find easily that all 
M a := (n a (x, t)) steady state = as long as all T a > 0, meaning the vacuum is absorbing for each color. As soon as 
some of the temperatures become negative, stationary solutions with M a > emerge, satisfying the corresponding 
equations g a (3Mp — —2r a . Thus, in general a multitude of hypersurfaces of first and second order transitions exist 
in the phase space spanned by the relevant temperature variables {r a } , which separate the phases where a specific 
color becomes extinct. Whenever Q2p QapMp + 2r a ) changes from a negative to a positive value, the order parameter 
M a undergoes a continuous or a discontinuous phase transition from an inactive absorbed state with M a = to an 
active state with M a > 0. All hypersurfaces of phase transitions meet in the multicritical point where all temperatures 
are zero. All homogeneous states are globally stable because the evolution of the total particle density of homogeneous 
states in time is given by d(J2 a M a )/dt = — J2 a ^aT a M a ~ | ^ Q ^ \ a g a pM a Mp. Thus all solutions of the mean field 
equations of motion are bounded to a finite region in the space of positive M a as long as the stability condition 
mentioned in the foregoing paragraph holds. 

In the following we focus on the effect of fluctuations on the scaling behavior of correlation and response functions 
in the vicinity of the multicritical point where the strongly relevant parameters {r a } are small. In order to apply 
field-theoretic methods and the renormalization group equation in conjunction with an e-expansion about the upper 
critical dimension |29|-|33| , i t is convenient to use the path-integral representation of the underlying stochastic processes 
n(x,t) = {n a (x,t)} |33| [36|. With the imaginary- valued response fields denoted by n(x, t) — {n Q (x, t)}, the generating 
functional of the connected response and correlation functions, the Greens functions, takes the form 



W[h,h] = In J V[n,n] exp(-j\n,n] + J d d x J dt(hn + hn) 



(4) 



The response fields n (x, t) correspond to the conjugated auxiliary variables of the operator formulation of statistical 



dynamics by Kawasaki |37| and Martin, Siggia, Rose P8| ). The dynamic functional j\n, n] and the functional 
measure I'm, nl, which in symbolic notation is proportional to Y[ x t {dn(x, i)dn(x, t)j , is understood to be defined 
using a prepoint (Ito) discretization with respect to time. The prepoint discretization leads to the causality rule 
8(t < 0) = in response functions which forbids response propagator loops in the diagrammatical perturbation 
expansion 



The Langevin equations (UN3) are recast as a dynamic functional 
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J = dtd d x ( n a (d t - A Q V 2 + R a (n)jn a - n a D af i(ji)np ) 

^ a a/3 ' 



= J dtd d x ^2 ^an a ^A Q 1 d t + r a - V 2 + ^ 9apnp ~ 9an a ^j ^ n Q . (5) 

In the second line we have neglected subleading terms. Correlation and response functions can now be expressed as 
functional averages of monomials of the n a and n a with weight exp(— J). A glance at Eq. ^ shows that the responses 
are defined with respect to additional local particle sources h a (x, t) > in the Langevin equations (Q). A rescaling 
of the fields n a — > l a n a , n a —> l~ x ri a leaves the functional J forminvariant but transforms the coupling constants 
9a/3 —> lfj9ap and g a — ► l~ l g a - Thus invariant coupling constants are given by f a p = g a pgp- A suitable rescaling is 
defined by g aa — g a . If we choose this normalization, we denote the rescaled fields by s a ~ n a and s a ~ n a . 

The scaling by a suitable mesoscopic length and time scale, fi^ 1 and (A/x 2 ) (with A Q ~ A) respectively, leads to 
s a ~ s a ~ /^ rf / 2 , / a /3 ~ /i e where e = 4 — d . Hence c? c = 4 is the upper critical dimension. It is now easy to show 
that all neglected possible subleading terms in the expansion of R a (n) and the noise correlation, as well as higher 
gradient terms in the Langevin equation (|l|) and non-Gaussian and non-Markovian noise correlations, have coupling 
constants with negative /i-dimensions near the upper critical dimension. Therefore, under the renormalization group 
flow, they are renormalized to zero and we can safely neglect them because we are interested in the leading universal 
critical behavior. These couplings can be reintroduced if one is interested in corrections to scaling. In a renormalized 
field theory context |2^,|3(J , the ^-dimensions of the coupling constants are equal to the so called naive or engineering 
dimensions. All coupling constants with positive naive dimensions (the relevant couplings with respect to the Gaussian 
fixed point as the starting point of the perturbation expansion) need a renormalization because the corresponding 
vertex functions are the only ones which develop primitive divergencies in perturbation theory. Thus the field theory 
based on the functional J (Eq. (|^)) is renormalizable and the calculated scaling properties are universal for the full 
class of MDP-processes. 

III. RENORMALIZATION AND ONE-LOOP CALCULATION 

Now we are in a position to develop the perturbation theory in the inactive phase with all the r a > 0. To begin 
with we separate the anharmonic "interaction" terms from the dynamic functional J, Eq. (H), and retain only the 
harmonic ones: 



Jo = y dtd d x (^Sg (d t + X a {T a - V 2 ))s Q - h a s a - h a s a 



(6) 



Here we have included the external sources h a and h a - The Gaussian path integral J T>[s, s] exp(— Jo) = exp((/i, Gh)) 
involves only non-negative fields s a (x, t) > as long as h a (x.,i) > 0. It yields the propagators for the Fourier 
transformed fields s(q, t) = J d d x s(x, t) exp(— iq ■ x) etc., as 

(s a (q,t)s (q!,t')) o = (2^) d ,5(q + q'X/3G Q (q,i-0 

G Q (q,t) = 0(t)exp(-A Q (r Q +( Z 2 )t) (7) 

where the Heaviside theta-function is defined with 6(t = 0) = following from the Ito-discretization of the path- 
integral and ensuring causality. 

Besides the propagators, the anharmonic coupling terms in J, Eq. (|^), define the elements of the graphical pertur- 
bation expansion. They are depicted in FIG. 1 where an arrow marks a J-leg and we draw diagrams with the arrows 
always directed to the left (ascending time ordering from right to left). 

q 

* =G a (q,t-t') 

a 
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FIG. 1. Elements of the graphical perturbation expansion 



FIG. 1 shows that the color of a sdeg on the left side of a vertex is not annihilated: at least one s-leg on the right 
side displays the same color. Thus, going from left to right, i.e. backward in time, through a diagram, colors have 
only sources and no sinks. This property is a consequence of the existence of absorbing states in the model. The 
perturbation expansion of a translationally invariant field theory can be analyzed by the calculation of the vertex 
functions T^ a y ^y({q,uj}) corresponding to the one-particle irreducible amputated diagrams. Here the sets {a} and 
{/3} denote the colors of the amputated outer s-legs and s-legs respectively. Going backward in time we conclude from 
the conservation property of the s-colors that the colors of the set {a} must appear as a subset of {(3}. Therefore the 
only nonzero two and three point vertex functions are r QjQ ,, r QQjCt , and T a ^ a p. Moreover, another property follows 
directly from color conservation: the vertex function ^{ a } ,{/?} does not depend on parameters A , r, and g with colors 
other than the ones of the vertex functions itself. Thus, T a ^ a , T aa a , and T a ^ aa are only functions of the parameters 
A Q , t q , and g aa = g a and are in particular independent from the interspecies couplings g a fj with a ^ (3. Therefore 
these vertex functions are the same as the corresponding functions of the well analyzed one-species Gribov process. 
To calculate them one can set the interspecies couplings g a p = 0. In this case the model shows rapidity reversal 
invariance s a (t) <-> —s a (—t) from which follows the equality T aa ^ a = —T a aa . 

Now we are ready to consider the renormalization of the model. It is known that the perturbation expansion of a 
field theory with a momentum cutoff A develops divergencies if A — * oo p9|j30[ |. If the model is renormalizable one 
can absorb all these "primitive divergencies" order by order in a loop expansion in a suitable reparametrization of 
the parameters. Absorbing the primitive divergencies regularizes the full model. The primitively divergent vertex 
functions have nonnegative naive [i -dimensions. Here, they are r QiQ ~ /i 2 and T aa _ a , T a ^ a p ~ \i e , logarithmic at the 
upper critical dimension, (FIG. 2). 
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FIG. 2. Primitively divergent vertex functions 

Taking into account the general properties of the vertex functions found in the foregoing paragraph, we see that 
the following renormalization scheme renders the theory finite 
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(8) 



Here G e = T (1 + e/2) j (47r) d ^ 2 is a convenient constant. Instead of a momentum-cutoff regularization, we use 
dimensional regularization and minimal renormalization in the following. From the discussion above we learn that 
the renormalization factors z\ a ' with i = (s. r, A) and := zi aa \ which are determined from r ajQ and T aa a , are 
already known from the one-species Gribov process and depend only on u a := u aa . Thus Z\ a) = Zi(u a ,e). The new 

renormalization factors Z^ 01 ^ — Zi nt ({u} , {A} , e) with a 7^ (3 stem from the interspecies couplings. They depend 
only on the couplings u a , up, u a p, up a , and the ratio X a /Xp. 

We will now explicitly calculate the renormalizations to one-loop order. The primitively divergent one-loop diagrams 
are shown in (FIG. 3). 
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FIG. 3. Primitively divergent one-loop diagrams 

Using dimensional regularization, we express the contribution of the self-energy diagram, FIG. 3(a), as a function 
of external momentum and frequency, q and lu: 



3(a) 



(A a #a) 2 



1 



G 

As 



p iuj + 2X a T a + A Q p 2 + A Q (p + q)^ 



ilu q 
2-e ~ r K + ^2 



(9) 



Here, we have retained only terms linear in u> and q 2 . These are the terms that display poles in e = 4 — d > 0. 

To determine the primitive divergencies of the vertex functions we can set external momenta and frequencies to 
zero and use equal temperatures r Q = r > as infrared regularisators. The contributions of the three diagrams 
FIG. 3(c,d,e) add to 



3(c) + 3(d) + 3(e) 



2 



2 



p 2A^(A Q + A^)(r + p2)2 
1 



+A 



G 



2A Q (A Q + A /3 )(r + p 2 ) 2 
1 



(2A Q ) 2 (r + p 2 ) 2 , 



A Q + A/3 



9aaga 



(10) 



From the zero-loop contributions and the results of our short calculation, Eqs. (|9|,|10|), we obtain the (as yet unrenor- 
malized) one- loop vertex functions to the desired order in u> and q 2 : 



G P 



4e 

+ X a T a ( 1 



1 - -jzg aa g a T a 6,2 ) + A Q q 2 ( l - -^g aa g a T~ 



■e/2 



G e 



2e(l-e/2) 



-e/2 



and 



r - \ « i 1 G £ f Kg a pgp + \pgi3 a g a \ e i 2 \ 

1 a,a/3 - A Q 5 Q/ 3 I 1 - — I — h g a a9a I T 7 I 



(11) 



(12) 



An explicit calculation of diagram (b) of FIG. 3 demonstrates that indeed T aa a — — T a ^ aa if g aa = g a . 

To absorb the e-poles in the renormalization Z-factors we note that the vertex functions are renormalizcd by the 
scheme Eq. (g) as 



r 



— ( 7S ai ) ■ ■ ■ 7X a ™)\ 

Qi -Q„ — \ ^ s J 



-1/2 



r, 



(13) 



G 



Using again the renormalization scheme Eq. (H), we find the renormalized vertex functions from Eqs. (0,121) as 



T aa = iu>[ zj°> - g(M/r Q ) E/2 J + A a q 2 ( - || ( M /r a ) £ / 2 



and 



+ ^-(^-5^y^ T -^) (14) 

r„, tf r w = G- V'A>„, (z<"" - Y e ( A °"t + ^ + "° + 2 "«) ( " /t) ' / ') ' (15) 



Therefore the vertex functions become finite by choosing 



= 1 + 1 ( XaU f + + ^ + 2^ 



(16) 



up to higher orders in the coupling constants u. As anticipated, for a = (3 we have found the well known renormal- 
ization factors of the Reggeon field theory. 



IV. RENORMALIZATION GROUP ANALYSIS AND ASYMPTOTIC SCALING 



Next we will explore the scaling properties of multicolored directed percolation. Scaling properties describe how 
physical quantities will transform under a change of length scales. At the end of chapter II we have introduced the 
arbitrary mesoscopic length scale \i. The freedom in the choice of fj,, keeping the unrenormalized bare parameters 

|f Q , A a , g a , g a /3^, and - in cutoff regularization - the momentum cutoff A fixed, can be used to derive in a routine 
fashion the renormalization group (RG) equation for the connected correlation and response functions, the Green 
functions 



N a +N a 



G ^ N H{ X ,t})=(i[[i[ S (^\t^) n ?( X f,f>) 



a 1 



(17) 



We denote /^-derivatives at fixed bare parameters by c^lp. From \i c^| G^ N ' N ^ = and the renormalization scheme 
Eq. (§), which leads to G<- N '^ = U a ( Z s a) ) (Na+ " a)/2 G {W '^ } , we then find the RG equations 



G {JV,iV} = 



with the renormalization group differential operator 



Here we have introduced the Gell- Mann-Low functions 

Ca=M^| lnA Q = 7 i a) -7i a3 , 

Pa? = fl d,\ Q u a p = (-£ + 7 (« + 7 | Q) + 7 f - 7 i Q/?) K/3 

with P a — (3 aa , and the Wilson functions 



(18) 



(19) 



(20) 



7; r) = M ®Aq ln Z » r) > i = s, A, r, n, r = a, a/3 



(21) 
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The RG equations (|l8|) can be solved in terms of a single flow parameter I using characteristics. Following this method, 
flowing parameters are defined by the characteristic equations 

l^jjU a /3(l) = (3 a f}{u{l)) , U a p(l) = U af3 , 

l^X a (l) = Ca(«(0)Aa(0 = C(u a (l))Xa(l) , A Q (1) = X a , 

dl 

l -J7 f <xO>) = Ka(u(l))f a (l) = K(u a (l))? a (l) , f Q (l) = T a , (22) 

dl 

and the RG equations (O) of the Green functions become 



d m-^ N a + N a . 



G^^>({x, t}, {f(l)}, {u(l)}, {A(Z)}, ln)=0. (23) 



Here, the functions 7 = 7i a \ C = Ccn an d independent of the color a and the interspecies couplings. 

The flow equations (g2| ) describe how the parameters transform if we change the momentum scale [i according to 
fi — ► p,(l) = Ifi. Being interested in the infrared (IR) behavior of the theory, we must study the limit I — ► 0. In 
general we expect that in this IR limit the coupling constants u a p(l) flow to a stable fixed point u a p.* according to 
the first set of Eq. (p2|). In particular, the intraspecies couplings u aa — u a then flow to a color independent fixed 
point u a (0) — because (3 a (u) = (3{u a ) , and this Gell-Mann-Low function (3 is equal to the corresponding function 
known from the one-species Gribov process. Thus, the fixed point value u* is independent from any coupling to other 
species. 

The solutions of the second and third set of the flow equations (^3|) are readily found in terms of the functions 
u a (l). In the IR limit I <C 1 they have the scaling form 

X a (l) = l z - 2 A{ a) X a , f a (l) = f-^A^Ta (24) 

where the Aj are nonuniversal amplitude factors. The scaling exponents 



are already known from directed percolation. From Eq. ( |23| ) we find the solution in the IR limit as 
G^>({x, t}, {r}, {„}, {A}, fi) = Yl {A M)(N a +N a )/2 

xG { ^({x, *}, {1 2 ~ 1 / v A t t}, K}, {l z - 2 A x X}, If,) (26) 
with AT = ^2 a N a , N = J2 a Na, and the DP anomalous field scaling exponent 

rj = 7(«„) . (27) 
Omitting the nonuniversal amplitude factors A[ and taking into account dimensional analysis in space and time 

G^>({x, t}, {r}, {u}, {A}, /i) = vP+WGW&ttv*, A/A}, {r/ M 2 }, {u}, {c}, 1) (28) 



where c a = X a /X, we get by combination of Eqs. (26 28j) the asymptotic scaling form of the Green functions 

G^>({x, t}, {r}) = l( N +^+n)/2 G {N,N}^ {t// i/- }) . (29) 

This has the important consequence that all scaling properties of the DP processes remain unaffected by the introduc- 
tion of many colors. Moreover, all intraspecies Green functions are completely independent of the coupling between 
the species, which follows from the absorbing state conditions for each color. 

However, the interspecies coupling constants u a p with a 7^ (3 determine the properties of the interspecies scaling 
functions. Therefore we will now consider the consequences of the flow equations for these parameters. For this 
purpose we need the Gell-Mann-Low functions (3 a f) explicitly. From the last of Eqs. ( |20| ) we know that each of these 
functions begins with the zero-loop term — eu Q/ 3, and the higher order terms are determined by the Wilson functions. 
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These functions, the logarithmic derivatives of the Z -factors, are given by 7 = /z c^Iq In Z = J2a8 flapd Ual3 InZ. In 

minimal renormalization the Z-factors have a pure Laurent expansion with respect to e: Z = l + Y^ /s + Y^ /e 2 -\ . 

Thus recursively in the loop expansion the Wilson functions also have a pure Laurent expansion and, because they 
are finite for e — > 0, this expansion reduces to the constant term, i.e. all e-poles have to cancel in the logarithmic 
derivation. Hence, we obtain the Wilson functions simply from the formula 7 = — YlaB u apd Ual3 Y lyl \ Now it is easy 



to get these functions from the one loop results of the Z-factors Eqs. ( |16[ ). We find to this order 



is ^ 1 ix 8 ' 2 

( Q( 3) _ 1 I ^ a U a f3 + \pUj3a 



from which one finds the Gell-Mann-Low functions Eq. ( p0| ) as 

/ 3u a hup \ a U af 3 + \pU/3 a \ 

0afi = {~ £+ — + — + 2(X a + X,) ) Ua0 ■ (31) 

The Gell-Mann-Low functions of the intraspecies couplings, the "diagonal" /3-functions, follow as (3 a — (— e + 
3u a /2)u a giving the stable fixed point values it Q * = ?i» = 2e/3. With the help of Eqs. (|20| , p5| , p7|) , this stable fixed 
point leads to the well known one loop order DP exponents z = 2 — e/12, v = 1/2 + e/16, and 7/ = — e/6 . Using the 
fixed point values M a * to obtain the Gell-Mann-Low functions of the interspecies couplings of a pair of colors a ^ f3, 
we get 

Pa/3 — l—r H ^tt — — • (32j 



3 2(A Q + A /3 ) 

In addition to the unstable decoupled fixed point values u a p* = up a * = 0, the equations f3 a p = f3p a = are solved by 
a fixed point line 

2A Q 2A^ 4e 

u a p* + t — — ^up m = — . (33) 



A Q + A,3 A Q + A/3 3 

This equation is the key result of this Section. Clearly, however, the one-loop calculation cannot give us any information 
whether the degeneracy of all the points on this line is fundamental to our model or is lifted by higher loop corrections. 
Thus, we must proceed to two-loop order, presented in the following chapter. We note, however, that there are two 
special points of unidirectional coupling on the fixed line: 

2(X a +X ) 

U aj 3* = , M/3a* = r-r £ and 

6Aj3 

2(A Q + A/3) 

Up a * = , lta/3* = — £ . (34) 

In the next chapter we will show that for colors with the same flavor, meaning A Q = A^j, these two points are indeed 
the only stable points on the line. Because M a /3* = or up a * — are fixed point values at any loop order, we conjecture 
that unidirectionality is generic for the asymptotic behavior of coupled DP processes, irrespective of whether colors 
belong to different flavors or not. 

V. TWO-LOOP RESULTS AND CROSSOVER TO UNIDIRECTIONALITY 

From the analysis in the foregoing chapters we know that the interaction of two colors does not depend on the 
existence of other ones. Thus in this chapter we consider a coupled model of two colors a — 1 and 2 with the same 
flavor, i.e. equal kinetic coefficients A a = A and intraspecies couplings g a — g. In contrast, the temperatures r a may 
be different. Thus, the (unrenormalized) dynamic functional is now given by 

J = J dtd d x\(^ 1 ^X- 1 d t + T 1 -\7 2 + ^(s 1 -J 1 )y i 

+s 2 (x~ 1 d t + t 2 - V 2 + |(s 2 - s 2 ))s2 + 2-(«i.9i2 + S292i)s lS2 ) . (35) 
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Note that in the case g\2 = 0, the dynamics of species 1 completely decouples from species 2 and vice versa. It follows 
that T a , a p — 0, if g a p = and a ^ (3. 

The detailed calculation of the two-loop contributions is presented in Appendices A to C. Adding Eq. (B6) to the 
zero- and one-loop parts of the selfenergy Eq. (|ll|), we get after renormalization, using Eq. (|l3| ) and the scheme Eq. 



Ar a Z t 



+A q 2 yz x - 

+iu> Z. 



2e(l-e/2) 



• 2 3 

1 + u 

1 e 16 



8e 

(n 2 /T a y/ 2 u 

4e 



/13 3 
1 + U fe-16 



1 + u 



4e 16 



C T 



c n 



(v 2 /T a y u 2 



{^/r a Yu 2 



O(s u 2 ) 



(36) 



where the constants Cj are given in Eq. (B7), and u — u a . Note that in comparison with Eq. (|l4]) now the one- loop 
terms acquire renormalizations to 0(u) in order to be consistent in 0(u 2 ) of the perturbation expansion. These 
renormalized one-loop terms are needed to compensate non primitive singular terms proportional to ln(/i 2 /r) arising 
now by the e-expansion of r QQ . Of course only primitive UV-divergencies, which means here r- independent ones, 
have to be regularized by renormalization to make the theory finite. Eliminating the £-poles from Eq. ( |36| ) by the 
Z-i actors, we find 



Z s = l 
Z\ = 1 
Z T = 1 



u u 2 {7 9. 4\ ^ ; 

1 3 + - In - + Otu 3 ) 

4e 32eVe 2 3/ V ' 



a 



u 



8e 128e V e 
u w 
2e + 2e\e 



{13 31 35 ,4 
■ + — In 



4 2 3 
+ 0(u 3 ) . 



+ 0(u 3 ) , 



(37) 



In the same way as for t he se lfenergy we find the other vertex functions. Renormalizing Eq. ( |l2] ) and adding the 
two- loop contribution Eq. ( C13 ) we get 



-p /^-l/2 e/2\ —1/2 | rr(aB) 17 -1/2 



s/2 



2e 



(> 



3li 3u U a p + Ufj a \ U a p + U/3 a 



2e 16 



2s 



16e 



+ 1 hr + 



4 1 

e + 2 



2 



- 2W _ 3U U af 3 + M/jq 

e 16 4e 



1 J U af 3Uf3 c 



(38) 



where Z u 



_ y(aa) 



eventually find 



Here the e-expansion also shows that the non primitive divergencies ~ ln(/i 2 /r) cancel. We 



Z^ (i) = 1 + (6u + u af} + U/Ja ) 



1 

lfe 



36 19 

T ~ T 



- - 7 )uU/3a + I 1 )Wq/3W/3. 



I 8 9 01 4 'l 



1 3, 4 

e~2 ln 3 



(<J + «jla) I +0(u 3 ) 



(39) 



and in particular 
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We are now in a position to calculate the renormalization group functions Eq. (|20| ) from Eqs. (|3^^,^) with the 
help of Eq. ( pl| ) . They are given by 

u (2 4\9w 2 ^. 
^ = -4 + U~ ln 3M + ° (w) ' 

u (17 , 4\ u 2 ^ / 3 . 

^-8 + U- ln 3jl28 +0(M) ' 

3m / 7 , 4\ 35m 2 ^. 3 , 
K = y - (T7; + ln o ) — + °( u ) . ( 41 ) 

and 



10 3 7 128 



A*/3 = ( S + U + J (w Q /3 + W/3 Q ) - iuQ/3W/? Q - ^ In - (u 2 a £ 



u 2 



13a . 



97 , 4\53m 2 /3 . 4\3uu a/J 5 s ,\ 

106 + ln 3 -64" ^ 4 - ln 3 ~ 32 ™*> + ° {u ) ) ^ " (42) 



Setting B — (3 aa — 0, we find the nontrivial stable fixed point value 

2e / /169 53 , 4\ \ 2n 

3 ( 1+ (288 + l44 ln 3j £ J + °^ (43) 

from which the scaling exponents Eqs. ( p5| , p7| ) of the Green functions Eq. (|2^) follow to second order of the e -expansion 
as 

e / f 25 161 , 4\ _ 2 .\ 

77 = — 1 + 1 In - £ + 0(e 2 ) , 

1 6y \288 144 3 J J 

„ £ ( , ( 67 59 , 4\ n . 2 ,\ 

z = 2 1+ H n - e + O £ , 

12 ^ V288 144 2,) K ' J ' 

The first two expansions have been known for a long time from Reggeon field theory |59[^0| where a different definition 
of the exponents is used. The expansion of the exponent v was presented by the author in (jj. The order parameter 
exponent 8, which enters the scaling law for the mean particle number in the active state < n >= M ~ | t - | ^ , follows 
from Eq. Q as 

/5 = i/ V = 1 -6( v 1 -^^ ln 3j £ + 0(£ " ) J ■ (45) 

We now turn to the interspecies coupling constants u Qj a with a ^ 8. The fixed point values as the solutions of the 
equation Bap = 0, where B a f) is the Gell-Mann-Low function Eq. (ji^), are of three different types: 

1. the decoupled fixed point U12* = U21* = 0, totally instable for u = u„; 

2. the two unidirectional coupled fixed points 1*12* = 0, U21* = 2u* + 0(e 3 ) and U21* = 0, U12* = 2u* + 0(£ 3 ); 

3. the symmetric fixed point M12* = U21* = u* + 0(e 3 ). 

To discuss the stability and the crossover between the last two types we try an Ansatz of the form u a p — w + ve a p 
with 612 = —£21 = 1- We already know from the one-loop result that u a p is driven by the renormalization flow to the 
fixed line w = it* with a crossover exponent cf> w = e/3 + 0(e 2 ). Setting therefore w — u in /3 a f3 we get 
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a p(u, v) = (3(u) + S af 30 v (u, v) , 

f3 v (u,v) = -Q - | In |) (u 2 -v 2 )v = -0.017 (u 2 v 2 )v . 



(46) 



For u — u* this equation shows that the symmetric fixed point v* = is instable in contrast to the stable unidirectional 
coupled fixed points = ±u*. The solution of the flow equation Idv/dl — v (u*,v) leads to the crossover 



v{lf = 



1+ (u 2 /v 2 - 



(47) 



with the very small crossover exponent 





3, 4\ 2 






G 


in - it* = 


g- 






4 3/ * 




3 3/ 



(48) 



A qualitative picture of the flow of the interspecies couplings in the plane u = u* under renormalization is shown in 
FIG. 4. 




FIG. 4. Flow of the interspecies couplings under renormalization 

In addition to the four fixed points D (decoupled), S (symmetric), U (unidirectional), the topology of the flow is 
determined by the symmetry line U12 — U21 which acts in the first quadrant as a separatrix between the regions of 
attraction of the two unidirectional fixed points. There exists another separatrix at the border of these regions, given 
to first order in e by U12 + U21 = 0, where the flow is driven to the hatched line U12 + U21 = — 2it*. On the left of 



this line we have the region of instability I, in which the condition ^2 a 



la^p > for all positive n a is violated. 



Therefore we conjecture that interspecies couplings with U12 + M21 < ultimately lead to first order transitions. In 
FIG. 4 the fixed point line of the one-loop calculation is also shown. This line is the support of the slow crossover to 
the unidirectional fixed points which therefore describe the ultimate critical behavior of the MDP universality class. 

We conjecture that in the case of colors with different flavors, \ a ^ A/3, all properties are smooth functions of the 
ratio \ a /\p as long as this ratio is sufficiently close to one. Thus generalizing to different flavors, the topology of 
the renormalization flow displayed in FIG. 4 should only be smoothly deformed but not destroyed. In particular, the 
unidirectional coupled fixed points U should be stable also for A Q ^ A/3. 
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VI. SYMMETRIES AND GENERAL FIXED POINT PROPERTIES 



Having found the detailed fixed point structure of the model Eq. ( pq ) of two colors with the same flavor in the 
two-loop approximation, we will now investigate which results are valid to all orders of perturbation theory. According 
to the considerations of the third chapter, one demonstrates easily that for a ^ ft the vertex function Y aa p = if 
9a/3 — 0. Thus the lines of unidircctionally coupled models in FIG. 4, namely u\2 = or U21 — , respectively, are 
invariant under the renormalization flow. Trivially, the decoupled fixed point D is the intersection point of these lines. 
For gi2 — 521 the dynamic functional J Eq. ( |35| ) possesses the symmetry s\ <-> S2, ~Sx <-> Si, and n «-> T2 from which 
we find that I\i2 = r2,2i- Thus, these two vertex functions need the same Z-factor: Z12 = Z21. It follows that the 
symmetry line in FIG. 4, u\2 = U21, is invariant under renormalization. 

Is there a condition that determines the crossover line? The answer is yes. We change to variables corresponding 
to the total and relative particle numbers, respectively, 



si + s 2 
si - s 2 



S=\( Sl +S2) , 

c=i(Vi-s 2 ). (49) 

Such linear transformations do not alter the measure of the functional integrals, and the dynamic functional changes, 
in the special case t\ = T2 = r, to 

J = J dtd d x A + (r - V 2 ) + (§ + ^±^i) s - §?) 



-c[ A- X a + (r - V 2 ^ ' 9 ^ 



s s 

; / 2 

; ) ! r - \ - ) + - sj c - |c _ s - gscc 
9 512 + 521 \ ^ + 912-921^,2 (50) 

We see that in the case gi2 + 521 = 2g the dynamics of the total particle number decouples from the dynamics of 
the relative one in the sense that all vertex functions containing s-, but no c -legs are zero if c-legs are attached. In 
particular Ts -cc — and Tg -ss — — T^^.s, which leads to the renormalized interspecies couplings with W12 + U21 = 2u and 
determines the crossover line to all loop-orders. Note that on the symmetry line gi2 = 521 holds, and the functional 
J is then invariant against c «-* — c, c«-> —c . 

The different fixed point values of the interspecies couplings are now fully determined as the intersection points of 
the several invariant lines (see FIG. 4): the unidirectional lines U12 = 0, U21 — 0, the symmetry line u\2 — U21, and 
the crossover line u\2 + U21 — 2w. Thus, we find to all orders for the symmetric fixed point S: U12* = M21* = it* and 
for the stable unidirectional fixed points U: uyi* = 0, 1*21* = 2u* and U21* — 0, 1*12* = 2u», respectively. These are of 
course the results that we have found explicitly in the two-loop approximation. 

The unidirectionally coupled model, which describes the generic scaling properties of the MDP processes, exhibits 
another symmetry. Using g\2 = and (721 = 2g', we write the dynamic functional of this model in the form 



Ju = J dtd d x X (V 1 ^ + (n 



V 2 ) + f( S i-si) )*i 



+S2^A 1 a t + (r 2 -V 2 ) + |(s2-S2)js2 + S2^-ff + 5's2jsij . (51) 

Here we have introduced a further harmonic unidirectional coupling ~ a, which corresponds to an additional linear 
term: dri2/dt = ■ • • + Aerrai in the Langevin equation for species 2. This term was first considered by Tauber et al. 
p7[ in their study of the nonequilibrium critical behavior in unidirectionally coupled DP processes. Such a term does 
not alter the general renormalizations as long as the corresponding composed field S2S1 is treated as a soft insertion. 
Of course, a is a relevant parameter like the temperatures T a and needs its own renormalization factor Z a determined 
by 



a = Z a Z- x a , (52) 
alized vertex function with an 

define the multicritical point. 



in such a way that the renormalized vertex function with an insertion Y 2 ^(sisx) ^ s ^ n ^ e ' settings a = r± = t 2 = 
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From Reggeon field theory one knows a transformation called rapidity reversal: s(t) <-> —s(—t), which is broken by 
the DP transition to an active state. Here we generalize it to read 

si(t) -«a(-t) , Si(*) -» -«a(-*) - , 

*a(*) -» . s 2 (t) -> -«!(-*) + * 2 (-t) • (53) 

Under this transformation the functional J u changes to 



J u -»■ J dtd d x A f«i ^A- 1 ^ + fa - V 2 ) + | (si - sr)^ Sl 



! A _1 9 4 + (n - V 2 ) + |(s 2 - s 2 ) ) « 2 



+s 2 ^-<r + ti - t 2 + g's^j si + (g - g') ^sis 2 si - s 2 si + s 2 s 1 s 2 j j ■ (54) 

We learn from this relation that in the case g = g' , J u gains a higher symmetry at the multicritical point, which is 
not destroyed by renormalization. It follows in this case that a renormalizes in the same way as the r Q , which implies 
Z a = Z T , and that the renormalizcd couplings are related by — u 2 i* = u* at the fixed point, as we already know 
from above. From these considerations follows that the crossover exponent, which is defined by the scaling invariants 
<j/t^ , is given by $ = 1. 



VII. THE THHG MODEL 

Recently Tauber et al. j27| (in the following abbreviated by THHG) introduced a general unidirectionally coupled 
DP process with species-independent diffusion coefficients, that was motivated by a study of Alon et al. [^5| on a 
nonequilibrium growth model for adsorption and desorption of particles which displays a roughening transition. 

The THHG-model reads in our dynamic functional language as 



Jthhg - J dtd d x\ ^ (V 1 ^ + n - V 



2 , 9_ t 
2 



Sl - Si Si 



+s 2 (^A 1 d t + T 2 - V 2 + |(s 2 -s 2 )js 2 

-82 (^r + bV 2 - y Sl - / 2 S 2 + /2S1 + y S 2 j sA 



(55) 



In contrast to THHG, we have introduced an additional cross diffusion term ~ s 2 V 2 si here that is indispensable for a 
complete renormalization of the general model. In physical terms, coarse graining will always produce cross diffusion 
in this coupled model. But coarse graining does more: it also produces a term proportional to the time derivative 
of the density of the first species in the Langevin equation of the second one - besides further irrelevant couplings. 
Accordingly, in the functional Jthhg a term proportional to s 2 9fSx arises. However, such a term can be eliminated 
by a suitable redefinition of the fields 

si — ► si - f3s 2 , si -> si , s 2 -> ? 2 , s 2 -> s 2 - /3si , (56) 

so that the harmonic parts with the time derivatives in the dynamic functional remain diagonal. Note that in the 



special case b = fx = f[ = f 2 = 0, f 2 = g' the functional Jthhg Eq. (55) is identical to J u , Eq. (|5l|), from which 
we know that it is fully renormalizable, and, in particular, in the case g' = g by the Z-factors Eqs P'ij^ol ). In the 
following we will demonstrate that the asymptotic properties of the THHG-model indeed belong to the universality 
class described by J u . 

It is well known that the infinitesimal generators of a continuous transformation of the fundamental fields, which 
leads to a forminvariance of the describing statistical functional, define redundant operators (composite fields). In 
the case that these redundant operators are relevant or marginal, they unnecessarily contaminate the renormalization 
group. Therefore they should be avoided from the outset. Here wc introduce a linear, homogeneous transforma- 
tion between the fields, which does not change the form of the functional Jthhg Eq. (p5[). Let us call it the 
a-transformation: 
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si — > si — as 2 , si — > si , s 2 — > s 2 , s 2 — > s 2 + asi . (57) 

Note the difference to the /^-transformation Eq. (56), which is exploited to eliminate a s 2 <9tSi-term from Jthhg- This 
a-transformation leads to new coupling constants that we denote with a bar: 

a = a + a (ti - r 2 ) , 

II = ft + 2af 2 + a (a - 1) g , f 2 = f 2 + ag , 

J[ = f' l -2af! 2 + a{a+l)g, J' 2 = f 2 -ag. (58) 
Now we renormalize the THHG-model by the scheme 

s\ -» si = Zs /2 (si + As 2 ) , si — ► si = Z s 1/2 si , 
s 2 -> s 2 = Z// 2 (s 2 + As x ) , s 2 -> ? 2 = Z// 2 s 2 , 
t q -> t q = Z a _1 Z t t q , g -> 5 = Zf^Z^Zgg , 
a -» a - Z^ 1 (Z CT (7 + (Fi + AZ r )n + (y 2 + AZ t )t 2 ) , 

/l -» A = Z^Z^ 1 (^i/i - ^ (! - ^) ^.9 - 2 ^^2/2 
/i - A = Z S - 1/2 Z^ (Z[f[ -A(l-A) Z g g 2AZ' 2 f 2 , 

h - A = Z.- 1/a ^A _1 (^/2 - ^) , 

ft -> A - ^ 1/2 ^ 1 fe/a - ^„s) • (59) 



where <? = (u/i £ /G e ) 1 / 2 , Z g = Zu^ 2 , and the Z^ with i = s, X,t,u are given by Eqs. ( |37j|40| ). Besides it, dimensionless 

coupling constants are defined as v a ^ — (G e / p e ) 1 / 2 f a \ The scheme Eq. (|9|) is chosen in such a way that the 
renormalized dynamic functional reads 



Jthhg = J dtd d x\ ^ (z.X^dt + Z T n - Z X V 2 + 

+s 2 (z^dt + Z t t 2 - Z A V 2 + ^ (s 2 - s 2 )^s 2 
+s 2 I 2AZ s \- 1 d t - Z a a - Y lTl - Y 2 t 2 - Z b 6V 2 



si — ai) si 



+M Sl + z 2 / 2 s 2 - Z^si - ^s 2 ) si j . (60) 

Note that the counter term 2AZ B serves to cancel primitive divergencies arising in the vertex function 9 w r~ . 
In dimensional regularization and minimal renormalization the counterterms A and Y a are given by series in 
beginning with simple poles ~ A^ /e and ~ Y^ /e, respectively. We have formerly shown that only the residua of 
these poles determine the renormalization group functions. 

It is now appropriate to define a-transformation invariant dimensionless coupling constants as 

Wi = y/u (vi + v 2 ) - v 2 , w[ = \/u (v[ + v' 2 ) — v' 2 , w 2 = Ju (v 2 + v' 2 ) . (61) 

The somewhat lengthy but simple calculation of all the one-loop renormalizations leads to the Gell-Mann-Low 
functions of the renormalization group equation. In part, in the case 6 = 0, they can be derived from the results of 
THHG [ p7| . If we set u to its fixed point value u* and define as usual j3 p = <9p/91n/^| , where p is any of the coupling 
constants, we obtain for the (5 -functions of the invariant couplings 

4(2wi + w 2 — u*)w 2 + 4:U r w[ + 16w*a — u* b(8w! + I0w 2 ) + 9(w*&) , 

&Pw' — 4:(2w[ +w 2 — u*)w 2 + 4w*wi + 16u»a — u*b(8w[ + 10w 2 ) + 9(u»&) 2 , 

4/3u, 2 = 8(w 2 — u*)w 2 + 4u*(u;i + w[) + 8u*a — 6u :t bw 2 + 3(w»6) 2 , 

8u*f3b — u*(w\ + w'i) + (w 2 — u*)w 2 + 16w*a — u*b(w 2 + w*) + (u*fe) 2 . (62) 
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The function a results from the additive renormalization A and mixes the renormalized fields in a correlation or 
response function under application of the renormalization group as 



V 4- — 



{si,si,s 2 ,s 2 } = {-as 2 ,0,0,-asi} 



Here 2? M is the renormalization group differential operator now given by 



T^n = Lidfj. + C^d\ + K T rd T + (K a a + Km + K 2 T 2 )d<j + 2J P P d p , 



(63) 



(64) 



with K(j — *y\ — 7 ff , k q = — y Q — a, and Eq. ( p3| ) acts on Green functions. The y a — 
additive renormalizations Y a . The function a — — ^2 p d p A^ is found to be 



a - 

The new renormalizations yield 

2/1 

.'72 



— f 2wi + mj + -l)w 2 - 4bw 2 + 3u6 2 



^p^yi 1 ' result from the 



(65) 





3b 




b 




v 2 v' 2 - 




fuv 2 - 




fuv' 2 




36 




b 




v 2 v' 2 - 




fuv' 2 — 


2" 


fuv 2 



3b 2 
4 

36 2 
4 



(66) 



In order to determine the fixed-point solutions of Eqs. (|62|), /3 P * = 0, we impose the condition a* = 6* = 0. This 
yields wi* = = 0, with w 2 * = (unstable) or w 2 * — u* (stable). It can easily be checked that these solutions 
are consistent with the full set of Eqs. (|6|) and (|65|). These are the solutions found by THHG J27j. Note that on the 
fixed point lines generated from the stable fixed point by the a-transformation a minimally coupled fixed point with 
v i* = v' u = v' 2lt = 0, w 2 * = y/iu is found. 

Now one has to prove stability of the fixed points of the full equations ( p^ ) without using the constraints a* = 6* = 0. 
A linearization about w\ = w{ = b = and either w 2 — or w 2 — shows that the flow of b is unstable for w 2 = 0, 
whereas it shows full stability of the fixed point for w 2 = u* . This vindicates the neglect of a, b and the corresponding 
counterterms A and Z^b in |27|| but only at the stable fixed point line, which is generated from the fixed point by the 
a-transformation. However without further knowledge this statement is only correct in the one-loop calculation and 
could be violated in higher loop orders. We will show that the stable fixed point is given by wi* = w' lif = 0, w 2 * = it* 
to all orders of the loop expansion. As a consequence, the fixed point values a*, 6*, yi*, y 2 * are zero, and j a * = 7r*- 
This leads to a crossover exponent $ = 1 where $ determines the scaling of ct/|tj|*. 

Indeed, we see from Eq. (61) that the stable one-loop order fixed point belongs, up to an a-transformation, to the 
dynamic functional Jthhg, Eq- (|o|), with coupling constants fx — f{ — f 2 = and f 2 = g, i.e. model J u , Ecl J55|), 
with the additional constraint g' = g. Above it was shown that this equality leads to rapidity reversal Eq. ( |53j ) as a 
higher symmetry. This higher symmetry is preserved under renormalization and, because J u is fully renormalizable, 
we have the result Wu = = 0, w 2 * = m* and $ = 1 to all loop orders. 

Computations based on the dynamic functional J7u, Eq. ([55]), are much easier to perform than calculations using 
the complete model Jthhg, Eq. (|60|). Thus it may be possible to find the equation of state for M 2 = g(s 2 ) to second 
order and check the assumptions made in |27j] on the reexponentiation of logarithms to yield the new order parameter 
exponent f3 2 of that paper (for a calculation of the equation of state for M\ = g{si) to two-loop order see [p~8|| ) . 

The model J u describes the coupled DP processes near the multicritical point t% = r 2 = a = 0. What is needed for 
a thorough calculation of (3 2 is a theory that comprises the limit a — + oo. Therefore our considerations here do not 
solve the problem addressed in |27| , namely the determination of the scaling exponent f3 2 that controls the scaling 
M 2 (x (jP 1 -! 32 (r lc — T\f' 2 where species 1 is in its active phase. THHG calculate (3 2 by reexponentiation of logarithms. 
(We have done a recalculation and find a slightly different value /3 2 = 1/2 — 13e/96 + 0(e 2 ). The difference arises 
from a subleading term resulting from the ominous peculiar diagram FIG. 9(c) in p7|.) However the approach of 
THHG relies on the assumption that simple reexponentiation is possible. To derive such scaling properties faithfully, 
one indeed has to solve the crossover problem a — > oo which (possibly!) induces a new scaling at infinity for the 
correlations of species 2. Some features of this crossover remind us of the crossover from special to ordinary behavior 
in the theory of surface transitions f43| , with a corresponding to the surface enhancement c and the species 1 and 2 
corresponding to the bulk and surface respectively. The crossover problem of interest here is thus as yet unsolved. 
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VIII. DISCUSSION 



In this paper, we have studied multicolored directed percolation processes (MDP). In particular, we have shown that 
the scaling behavior of these coupled DP processes near their absorbing state transition is determined by the same 
critical exponents as known from simple one-species DP and therefore is independent from the number of colors. The 
characteristic asymptotic feature of MDP shows up in the asymptotic unidirectional coupling of each pair of colors. 
A special result of our analysis is the very slow crossover to this asymptotic unidirectionality which may be seen in 
computer simulations. The unidirectional behavior of the couplings of an interacting population is summarized in the 
following graphical picture. Consider a graph where each node represents one color. The stable fixed points are then 
represented by the so-called tournaments, that are the complete graphs with directed edges |Q. The directed edge 
from color a to color (3 stands for the influence of a on (3 in the respective equations of motion. In particular, for 
a population of three species there exist two different tournaments (up to permutation of the colors) which we call 
"cyclic" and "hierarchic", FIG. 5. 





2 1' — * — »2 

(a) (b) 

FIG. 5. Cyclic and hierarchic tournament of three species 

One gets a graphical picture for the not fully stable fixed points either by deleting the directionality of the edges 
between the nodes for symmetric couplings of the corresponding colors, or by completely deleting the edges for the 
uncoupled pairs (uncomplete tournaments). 

To get a simple qualitative impression of the behavior of the dynamic system of a population corresponding to a 
tournament of colors with the same flavor, we consider, in the active region with allr Q := — r a /2g > 0, a renormalized 
mean-field theory of the equations (|l],|^) for spatially homogeneous densities n a and set g a p = g (1 + e a p) with e aa = 0. 
We redefine the time scale t — > 2t/Xg and get 



d t n a = 







(67) 



Then the decision, which of all the e a p = —Ep a equal ±1, defines the tournament. The edge between a pair of species 
is directed from [3 to a if e a p = —Ep a = +1 and vice versa. The directions of the edges therefore represent the 
unidirectional "pressure" on the reproduction rate resulting from one color to another. Despite the simplicity of the 
Eqs. (^Vthey can generate a complex dynamic behavior (see e.g. (2^j24|]). First, let us consider the stationary states 



of Eqs. (pTft for a population consisting of three species. In the ternary phase diagram spanned by the positive rates 
?"i)7*2>7"3 m the subspace ^ Q r a = r = const., one finds different regions with one, two ore all three species alive, 
FIG. 6. 




/ 1 






\ 1+2 \ 



(a) (b) 

FIG. 6. Phase diagrams of the cyclic and hierarchic tournament 
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These regions are bounded by critical lines with absorbing state transitions where some colors become extinct. 
The dynamic behavior of the hierarchic tournament FIG. 5(b) is relatively simple: from each nonequilibrium initial 

state , the system relaxes to the stationary state which is a stable node. In the case of the cyclic tournament 

FIG. 5(a), for rates {r a } such that we have a three species stationary state, one also finds regions for which the ultimate 
relaxation behavior is characterized by attracting nodes. But for rates which lead to stationary states belonging to 
the crosshatched area in the ternary diagram FIG. 7 spanned by the we find stable spirals (damped cyclic 

relaxation) as the ultimate relaxation behavior (qualitatively pictured by the trajectory in FIG. 7). 




FIG. 7. Dynamic behavior of the cyclic tournament 

For stationary points near the middle of FIG. 7 the damping is very small and cyclic behavior dominates the 
dynamics. Especially in the middle of the diagram, i.e. for equal rates n = r2 = ra = r/3, the motion is not damped 
anymore. In this case, the dynamic system Eq. j67|) is known as a special form of the May-Leonard model Ef| that has 
been extensively studied in mathematical biology p3| , [l6| -ff8[ . Finally, after a relaxation in the plane n — ^ Q n a = r, 
there exists another constant of motion m = Y[ a and the dynamic behavior is characterized by limit cycles around 

the neutral stationary point 

Summarily we have shown that also in stochastic multispecies models of populations that evolve near the extinction 
threshold of all colors and therefore have many absorbing states, the critical properties at the multicritical point and 
at all continuous transitions are governed by the well known Gribov process (Reggeon field theory) exponents (for a 
previous simulational result on a two-species system which seems to agree with our findings see |49(| ) . In other regions 
of the phase diagram of course more complicated critical behavior may arise such as multicritical points with different 
scaling exponents [^o| . The models considered here have many absorbing states (each combination of colors may go 
extinct irrespective of the other ones), and therefore are different from models which were considered by Grinstein 
et al. []l4| where it was shown that multispecies systems with one absorbing state belong to the Gribov universality 
class. In addition we have shown that the universal properties of interspecies correlations and the phase diagram 
are determined by totally asymmetric fixed point values of the renormalized interspecies coupling constants. This 
eventually leads to a system working cooperatively. It is interesting that the asymmetry between the species seems to 
be the condition for this cooperation near extinction. The model considered here is a simple but universal model of 
such a cooperative society and should therefore have many applications in all fields of natural and even social science. 

Competition and extinction is of course a subject much considered in theoretical biology. The main differences 
between the present work and the topics covered e.g. in the monograph of Hofbauer and Sigmund are the more 
realistic local description of the interactions between the species, their diffusional motion, and the inclusion of local 
fluctuations. Thus here the equations of motion are local stochastic partial differential equations. Coarse graining 
and renormalization lead to an universal macroscopic picture of cooperativity near the critical states of extinction. 
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APPENDIX A: TWO-LOOP INTEGRALS 



In the calculation we will encounter momentum integrals of the type 



G. 



1 



qi i (<7i 2 + t) (q$ + r) (V + g 2 2 + (<& + <z 2 ) 2 + 3r) 



(Al) 



where G e = T(l + e/2)/(47r) rf / 2 , e = 4 - d, and / . . . = (27r)~ d / d d q . . . . They can be derived from two "mother" 
integrals 



M^(a,b;c) = G e 



M (2) (a;c) = G £ 



?1 + a )[<li + b)lq? + qi + (qi + q2? + C 



91,92 



(?l + a) (<Z 2 + 9 2 2 + ill + 92) 2 + c) 



(A2) 



by taking derivatives with respect to the parameters a, 6, c. Discarding nonsingular terms, we find in dimensional 
regularization 



M^(a,b; c) 



(A3) 



These formulas yield the singular parts of the integrals (Al) as 



r (SP) 

'll;l 

r(SP) 
'l2;l 



2 + 3e 



r (SP) 

! 11;2 



2e 2 



3eln- 



t (sp) J_, 2 _ £ + £ln 



£ 3 

r(SP) _ 

'10:3 — 



T-(SP) 

*30;1 



12e 



8e 2 



r(SP) 
'20:1 



2c 



In- 



(A4) 



APPENDIX B: TWO-LOOP SELFENERGY DIAGRAMS 





(a) (b) 

FIG. 8. Two-loop selfenergy diagrams 
In FIG. 8 the two- loop selfenergy diagrams are drawn. Diagram FIG. (8a) leads to 

3 



4(a) 



(A.9) 4 



J J fl dUe-^ G(q/2 + qi , ii)G(q/2 - qi, ti - h)G(q/2 - qi , f 3 ) 



qi,q2 



xG(q/4 + q 2 , ti - i 3 )G(q/4 - q x - q 2 ,ii - t 3 ) 



(Bl) 
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where G(q, t) is the propagator (Eq. (Q)) and we always set r Q = r > as an IR regulator. Noting that e zwtl G(ii) = 
(e"^ (tl ~ t2) G(ii - i 2 ))( e - iu ( t *-t 3 )(}fi 2 _ t 3 )) (c" iw * 3 G(t 3 )) factorizes in the integral, the time integrations over the 
intervals (ti — t 2 ), (t 2 — is), £3 are easily performed. The expansion in iu> and q 2 to linear order eventually yields 



4(a) = ^G E 2 r- e fr/ 20;1 - y (/ 30; i + 



{laB;l + T^20;2 + 7Tj(^20;3 — ^10;3) ) ) (B2) 



2 V ou,i 4 ' 2d 



where the Iki- lT n denote the integrals defined in Eq. (Al). Extracting the singular parts using Eq. (A4), we obtain 

4(fl) =32^ T r + (e +5 + ln 3)2A + U + 12 + 2 ln 3)Tj- (B3) 
In the same way we calculate the second selfenergy diagram FIG. 8(b): 

3 

4(o) = (A 5 ) 4 f f 11,1^-^0(^2 + q 1 ,t 1 -t 2 )G(q/2-q 1 ,t 1 -t s ) 

i=l 

xG(qi - q 2 , ta - t 3 )G(q/2 + q 2 , t a )G(q/2 - q 2 , t 3 ) 
- y-G 2 T~ £ (tI U;1 - y (j 12;1 + 7 11;2 ) 

-y (il2jl + /ll;2 + ^( 7 H;3 + 27l0j3 - Zll; 2 ))) (B4) 

up to higher orders in iu and g 2 . Extracting the singular parts again, we find 

Summing up, we finally get from 4(a) and (46) the two-loop contribution of the selfenergy as 

ri 2-ioo P ) = hE^l T -e ( CtT + Cj^ + CfA , (B6) 



where 



1 9e ^ 7/ 3e 9e, 4\ „ 13/ 19e 35e , 4\ 

G T = — I , C^ = — H In - . G„2 = H In - . (B7) 

2 16' 32 V 14 14 3/ 9 128 V 52 26 3/ V ; 



APPENDIX C: TWO-LOOP VERTEX DIAGRAMS 



(a) (b) (c) 

(d) (e) (f) (g) 

(h) (i) U) (k) 
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FIG. 9. Two-loop vertex diagrams; the numbers show the other time orderings 

FIG. 9 presents the eleven two-loop vertex diagrams. They can be calculated by the same method as the selfenergy 
diagrams. Here the external frequencies and momenta can be set to zero because the expansion in these variables 
does not lead to primitive divergencies. As an example we show the calculation of the diagram FIG. 9(i) explicitly. 

The two possible different time orderings of the vertices are shown in FIG. 10. 





FIG. 10. Two different time orderings of a vertex diagram 

The symmetry factor of the diagrams is one, but one also has to add the two diagrams arising from the interchange of 
the color indices f3 and 7. After the factorization of the propagators, the integration over the time intervals, indicated 
by the broken lines, is trivial und leads to the expression (with the abbreviations /Cj = A(r + q 2 ), = qi + q2) 



5(i) 



A* 

' 8 



1 
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Thus we obtain the contribution of the diagram FIG. 9(i) to the vertex function T a ^ a p: 
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In the same manner we calculate the other diagrams of FIG. 9 and find 
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Adding up all the two-loop contributions r^jL, • • • r^^g and extracting the singular parts Eq. (A4), we finally get 
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